Instability and entanglement of the ground state of the Dicke model 
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Using tools of quantum information theory we show that the ground state of the Dicke model 
exhibits an infinite sequence of instabilities (quantum-phase-like transitions) . These transitions are 
characterized by abrupt changes of the bi-partite entanglement between àtoms at critical vàlues Kj of 
the atom-field coupling parameter k and are accompanied by discontinuities of the first derivative of 
the energy of the ground state. We show that in a weak-coupling limit (ki < k < «2) the Coffman- 
Kundu-Wootters (CKW) inequalities are saturated which proves that for these vàlues of the coupling 
no intrinsic multipartite entanglement (neither among the àtoms nor between the àtoms and the 
field) is generated by the atom-field interaction. We analyze also the atom-field entanglement and we 
show that in the strong-coupling limit the field is entangled with the àtoms so that the von Neumann 
entropy of the atòmic sample (that serves as a measure of the atom-field entanglement) takes the 
value Sa = 5 ^ n (N + 1)- The entangling interaction with àtoms leads to a highly sub-Poissonian 
photon statistics of the field mode. 



Interactions between quantum objects lead to correla- 
tions that may have no classical analogue. These purely 
quantum correlations, known as entanglement, play fun- 
damental role in modern physics and have already found 
their applications in quantum information processing and 
communication [1]. A degree of quantum entanglement 
depends on the physical nature of interacting objects and 
on the character of their mutual coupling. There is an- 
other physical phenomenon that shares many of the fea- 
tures with the entanglement - quantum phase transitions 
(QPT) [2]. These transitions occur at zero temperature 
and they are induced by the change of a coupling param- 
eter. Recently several authors have reported on a close 
relation between quantum critical phenomena and ap- 
pearance of entanglement. In particular, scaling of entan- 
glement close to a quantum phase transition in an Ising 
model has been reported [3]. In this Letter we report 
on an infinite sequence of instabilities (i.e., a sequence of 
quantum-phase-like transitions) and its connection to bi- 
partite atòmic entanglement in the so-called Dicke model 
(DM) [4-8]. The QPT in the DM was first rigorously de- 
scribed by Hepp and Lieb [9] (see also Ref. [10]) while the 
instability of the ground state has been reported by Nar- 
ducci et al. [6]. In the paper by Hepp and Lieb [9] only 
a thermodynamics limit of large-7V and weak-coupling k 
has been studied. It has been shown that in this limit 
at k — u>/ y/N the ground state of the model exhibits 
the second-order QPT. QPT in the DM have been stud- 
ied by Hillery and Mlodinow [11] and Rzazewski et al. 
[12] in the limit of large N and weak coupling. QPT 
and its relation to a chaos have been studied by Emary 
and Brandes [13]. Schneider and Milburn [14] (see also 
Refs. [15,16]) studied appearance of entanglement and 
its connection to quantum phase transitions in the DM 
with just two àtoms. The connection between the atom- 
field entanglement and the QPT in the thermodynamical 
limit of the DM has recently been studied by Lambert 



et al. [17]. However, in all these studies the sequence 
of ground-state instabilities (quantum-phase-like transi- 
tions) and the behavior of bi-partite atòmic concurrence 
have been completely overlooked. In this Letter we re- 
port on this remarkable feature of quantum entanglement 
in the Dicke model both in the weak-coupling as well as 
strong coupling regimes. 

Let us assume a model of N identical two-level àtoms 
(spin-1/2 partides), at positions ri, . . . , r/v coupled to a 
single mode electromagnètic field via the electric-dipole 
interaction originally proposed by Dicke [4] (see also 
Refs. [5-8]). In the rotating-wave approximation the 
model Hamiltonian reads 

H = A cr| + Tilúfo* a 



3=1 

N 
3 = 1 



(1) 



where a and arc the field annihilation and creation 
operators, respectively. In what follows we will assume 
that the atòmic frequency u>a is resonant with the field 
frequency up and h = 1. All the àtoms are assumed to 
be coupled to the field with the same coupling k and 
the position dependent phase factors can be included 
into defmition of the Pauli matrices <j via local uni- 
tary transformations (see e.g. Ref. [11]). As shown by 
Tavis and Cummings [5] one of the integrals of motion 
of the DM corresponding to the total excitation numbcr 
P = a + a j a J partitions the total Hilbert space 

of the DM into a direct sum of subspaces labelled by 
different excitations numbers. 

When the number of excitations p is smaller than 
the number of àtoms N, then the corresponding sub- 
space of the Hilbert space is spanned by p + 1 vec- 



tors \{e s ,g 



* n ®{N-s) 



})a\p — s)f with s — 0, . . . ,p where 
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|{e s , g®( N ~ s ^})A is a completely symmetric state of s 
àtoms in the upper level |e) and (N — s) àtoms in the 
lower level \g), while \p— s)f describes a Fock state with 
(p — s) photons. The DM Hamiltonian (1) in this sub- 
space is represented by a 3-diagonal (p + 1) x (p + 1) 

matrix the eigenvalues (k) of which are linear func- 

tions of the coupling k, i.e. E^ p \k) — K^k + Eq P \ 

where Eq P ^ = (jp — N/2)uj and are parameters to be 
determined from the eigenvalue problem. The lowest en- 
ergy eigenstate in the givcn subspace of the Hilbert space 
is determined by the smallcst (negative) K^J n . In what 
follows we will omit the subscript min and when we use 
EW(k) and we refer to the energy of the ground 

state for a given value of k. 

The character of the ground state of the DM (1) with 
the smallest eigenenergy, depends on the coupling n. At 
the critical vàlues of the coupling kj the number of excita- 
tions p in the ground state of the system, the bi-partite 
concurrence C between àtoms as well as the von Neu- 
mann entropy of the field mode do exhibit discontinuities. 
These effects can be clearly seen from Table 1, and Figs. 1 
and 2. Specifically, let us consider cigenenergies of the 
ground state with and 1 excitations. Thcir cxplicit ex- 
pressions are EÍ°) = -f Tiw and = ^-cj-ny/N. It 
is clear that for k/uj < 1/y/W the energy EÍ°) is smaller 
than Ei 1 ), so the state \g® N ) with zero excitations is 
the ground state of the DM. But for k/uj > Í/VN this 
state is not a ground state of the DM anymore since 
Ei 1 ) < Ei°). That is n x = lj/Vn is the critical point 
at which number of excitations in the ground state of 
the DM is abruptly changed. Analyzing the inequali- 
ties between eigenenergics of lowest energy eigenstates 
in subspaccs of the Hilbert space with different excita- 
tions numbers we can determine critical vàlues of the 
coupling parameter k. The first few critical vàlues can 
be expressed in a simple analytical form 

«i/w - N- 1 ' 2 ; (2) 
k 2 /uj= [V4JV-2- Vn]' 1 ; 

k 3 /u = [\J 5(N - 1) + y/(W - 5) 2 +8N- V^N - 2}- 1 . 

In a thermodynamical limit of large N (while p <C N) all 
critical points "merge" together, i.e. Kj ~ lú/V~N (where 
j = 1, . . . ,p) and we recover the result of Hepp and Lieb 
[9]. 

On the other hand in the strong-coupling regimc (p > 
N) when the Hamiltonian is represented by the (N+l) x 
(N+ 1) matrix (the size of which does not depend on the 
particular value of p) the energy of the ground state of the 
DM again lincarly depends on k and for a givcn p 3> N it 
can be approximated as EÍ p ) = —N^/pn + peu. The crit- 
ical vàlues of n in the limit p ^> N can be approximated 
as k p ~ üj[N(v / p+ 1 - ^/p)] -1 . 

In our approach we consider distinguishable àtoms (lo- 
cated at positions f\, . . . , rjy) that interact with a single 



mode field. This field mediates quantum correlations be- 
tween the àtoms. In what follows we will study bi-partite 
atòmic entanglement in the ground state of the Dicke sys- 
tem that for a given k is represented by a superposition 
= ELo A s\{^ s ;g^ N ~ s) })A\p- s) F . The com- 
posite states \{e® s ; g®( N ~ s )}) of this superposition ex- 
hibit entanglement. For instance the state |{e; g®( N ~^}) 
is the state (the so-called entangled web) with maxi- 
mal bi-partite entanglement between an arbitrary pair 
of àtoms (see Koashi et al. [18]). 

In order to make our discussion quantitative we re- 
cali the definition of the concurrence [19,20] between 
a pair of qubits (two-level àtoms). Let g is a den- 
sity matrix of the pair of two-level àtoms expressed in 
the basis {\gg}, |<?e), \eg), |ee)}. Let g is a matrix de- 
fined as g = (a v ® a y )g T (cr v ® a y ), where a v is a 

Pauli matrix ^ ? ^ ^ , while T indicates a transposi- 

tion. Thcn the concurrence of a bi-partite system reads 
C = max{Ai — À 2 — A 3 — A 4 , 0}, where Ai are the square 
roots of the eigenvalues of the matrix gg in descending 
order. The vàlues of concurrence range from zero (for 
separable states) , to one (for maximally entangled states 
of two qubits). The bipartite entanglement of the entan- 
gled web state |{e; g®^ -1 )}) measured in concurrence 
takes the value C = 2/N. As follows from our investi- 
gations (see Table I and Fig. 1) for K\ < k < n 2 the 
entanglement between an arbitrary pair of àtoms corre- 
sponds to the concurrence CÍ 1 ) — 1/N. Increasing the 
coupling the value of the concurrence is decreasing but it 
changes discontinuously at critical vàlues Kj (see Fig. 1). 
Nevertheless, the total amount of the bipartite entan- 
glement in the atòmic sample measured in terms of the 
function ta is non-zero for arbitrarily large number of 
àtoms and arbitrarily large coupling (see Fig. 2). 

The bi-partite entanglement between àtoms is medi- 
ated via a single-mode electromagnètic field that inter- 
acts with the àtoms in the dipole and the rotating-wave 
approximations. In fact, this interaction leads also to an 
entanglement of the field mode to the atòmic system. In 
order to measure the degree of entanglement between the 
field and the àtoms we utilize the von Neumann entropy 
of the atòmic sample Sa = — ^[pa Pa], where pa is 
the density operator of the àtoms in the ground state 
\EÍ p )). It follows from the Araki-Lieb theorem that for 
a pure atom-field state \'EÍ p ))af the field entropy Sf is 
equal to Sa- Larger the entropy is the larger the degree 
of entanglement. For p < N the upper bound on the 
entropy Sa is given by the number of excitations p, i.e. 

s (rnax) = ^ + ^ for p > ^ thfi 

upppcr bound 

is given by the total number of àtoms in the system, i.e. 

a given number of àtoms the 
atòmic (field) entropy depends on the number of excita- 
tions in ground state of the DM. For very weak coupling, 
i.e. when > k > k,\ the àtoms are totally disentangled 
from the field. In the region of vàlues K\ < k < k 2 when 
the ground state is characterized by one excitation, the 
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field mode occupies a Hilbert space spanned by two Fock 
states \0}f and |1)f- In this case, the entropy is maximal 
and equals to ln 2. With the increase of the coupling n 
the entropy of the field is always smaller than ln(p +1). 
In fact, for large N and for p large ( p > N) we find that 
Sa = \ ln(iV + 1), i.e. the entropy is equal to om half 
of its maximally possible valuc. That is the field is not 
maximally entangled with the atòmic samplc. The rea- 
son for this behavior is that the entanglemcnt cannot be 
shared frccly. To illuminate this issue let us consider the 
range of k such that p = 1. In this case the field mode 
can be effectively represented as a qubit. Consequcntly, 
the systcm of N two-level àtoms and the field mode can 
be represented as a set of (N + 1) qubits. Recently Coff- 
man et al. [21] have conjectured that for pure states the 
sum of square of concurrències between a qubit j and 
any other qubit k is smaller than or equal to the tangle 
between the given qubit j and the rest of the system, i.e. 
the CKW inequalities read 



N 

E 



r 2 



< T, 



3,3 



4detpj 



(3) 



where the sum in the left-hand-side is taken over all 
qubits except the qubit j, while r. j denotes the tangle 
between the qubit j and the rest of the system (denoted 
as j). If wc assume that the qubit j represents the field 
mode, then we can find that the tangle between the field 
and the system of àtoms reads tfa = 4dctp^ = 1, while 
the concurrències between the field and each of the àtoms 
is CfAj = 1/VN. From here it directly follows that the 
ground state of the DM with p = 1 saturates the CKW 
inequalities, which proves that the atom-field interaction 
as described by the Hamiltonian (1) with small coupling 
(«4 < k < k 2 ) docs induce only bi-partite entanglement 
and does not result in intrinsic multipartite quantum cor- 
relations. For the moment it is csscntially impossible to 
generalize this result for other vàlues of p since no mea- 
sures of entanglemcnt between a qudit (field mode with 
p > 1 excitations) and a set of qubits (àtoms) are avail- 
able and no corresponding gcncralization of the CKW 
inequalities is known. Neverthcless, it is interesting to 
study the strong-coupling limit whcn p 3> N from a dif- 
ferent perspective. In this case the ground state of DM 
reads \E^) AF = £f=o A s \{e®°;g®í N -ï}) A \p-s) F . For 
N <C p the amplitudes A s can be approximatcd as A s ~ 



N 



-N 



We already know that in this 



case the field mode is highly entangled with the àtoms 
that is reflected by the field entropy Sf — \ ln(N + 1). 
Quantum correlations that are establishcd in the atom- 
field system lead to highly non-trivial photon statistics 
of the field mode. Specifically, from above it follows that 
for p ^> N the mean excitation of the atòmic sample is 
equal to N/2. Correspondingly, the mean photon num- 
ber of the field mode is equal to n = p — N/2 with the 
dispersion of the photon number distribution equal to 



A = \i N /2. This means that in the strong-coupling limit 
entanglement between àtoms and the field leads to highly 
sub-Poissonian photon statistics of the field mode in the 
ground state of the DM. This is reflected by the Mandel 
parameter of the field mode, defined as Q = (A — n)/n, 
that is close to -1 for p 3> N. 

In summary, using tools of quantum information 
theory we have shown that the ground state of the 
Dicke model exhibits an infinitc sequence of instabilities 
(quantum-phasc-like transitions) . These transitions are 
characterized by abrupt changes of the bi-partite entan- 
glement between àtoms at critical vàlues Kj of the atom- 
field coupling parameter k. 
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TABLE I. The number of excitations p in the ground state of the DM with N àtoms changes abruptly at critical vàlues 
Kj. The larger the coupling is the larger the number of excitations in the ground state. The energy of the ground state E 
is continuously decreasing with the increase of the coupling re. For different regions of re characterized by different excitation 
numbers the ground-state energy is presented in Table 1 and in Fig. 1. At critical vàlues Kj the first derivatives of the energy 
OE/Ok exhibit discontinuities while the bi-partite concurrence reflecting entanglement between pairs of àtoms is changed 
abruptly (discontinuously) . 




FIG. 1. The energy E (measured in units of TuS) of the 
ground state of the DM with iV = 12 àtoms as a function of 
the scaled coupling parameter re/w. We see a sequence of in- 
stabilities (quantum-phase-like transitions) , that are marked 
by vertical lines corresponding to different vàlues of the cou- 
pling Kj . The first interval of vàlues of the coupling parameter 
< re < rei is characterized by the excitation number equal 
to zero. The second interval rei < re < K2 is characterized by 
one excitation, etc. The first derivative of the ground-state 
energy (black line) is not a continuous function of re (see the 
inset that illustrates the behavior of the energy of the ground 
state around the critical point Kj). In the figure we plot also a 
bi-partite concurrence as a function of re. The concurrence ex- 
hibits abrupt changes at "critical" vàlues of Kj . The maximal 
value of bi-partite concurrence in the Dicke system appears 
for rei < re < K2 when the ground state of the system contains 
just one excitation (p — 1). Then, the value of the concur- 
rence decrease with the increase of p. Nevertheless, even for 
p S> 12 the concurrence is non-zero. 



FIG. 2. The bi-partite concurrence of the ground state of 
the DM as a function of number of àtoms iV and the number 
of excitations p. Since the bi-partite concurrence decreases 
as 1/N we plot the "total atòmic bi-partite entanglement" 
ta = C 2 N(N - l)/2, where the factor N(N - l)/2 corre- 
sponds to the number of possible pairs in the atòmic sample. 
We see that ta as a function of N takes the largest value for 
p = 1 and in the large iV limit tends to the value 1/2. We see 
that even in the strong coupling case when p — N arbitrary 
pair of àtoms in the sample is entangled. This is reflected by 
non-zero value of ta (see the diagonal line in the (p, N) plane 
of the figure). In this case, for iV large enough the field ex- 
hibits sub-Poissonian photon statistics with the mean-photon 
number n = 2N/3 and the dispersion A = y / 5iV/26. We 
conclude that larger the number of excitations p smaller is 
the amount of bi-partite entanglement between àtoms. In 
the strong-coupling limit a strong atom-field entanglement is 
established that is quantified by the von-Neumann entropy 
S A = S F = im(AT + l). 
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